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ABSTRACT 


The  purpose  of  this  thesis  is  to  obtain  the  approximate  power  of 
a  test  for  serial  independence  in  a  stationary  linear  Markov  process , 

The  approximate  power  is  obtained  for  a  test  based  on  a  sample  estimate 
of  the  serial  correlation  coefficient  in  a  stationary  linear  Markov 
process  with  unknown  mean.  Chapter  I  introduces  the  problem  and  Chapter 
II  reviews  Patrick's  derivation  of  the  approximate  distribution  function 
of  a  sample  serial  correlation  coefficient  in  a  circular  universe.  Then, 
in  Chapter  III,  following  Patrick's  approach,  the  approximate  distribution 
function  of  a  sample  serial  correlation  coefficient  for  the  stationary 
case  is  derived.  This  is  used  in  Chapter  IV  to  obtain  the  approximate 
power  of  the  test  for  serial  independence. 
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CHAPTER  I 


INTRODUCTION 

The  recent  interest  in  the  problem  of  testing  for  serial  correlation 
in  linear  Markov  processes  seems  to  have  been  initiated  by  a  paper  of 
R.  L,  ANDERSON  (1942)  .  In  this  paper,  using  the  circular  definition  of  the 
serial  correlation  coefficient  as  introduced  by  Hotelling,  the  author  derived 
a  useful  test  of  significance  and  gave  a  table  of  significance  points  for 
the  lag  -  L  serial  correlation  coefficient.  A  summary  of  some  of  the  many 
papers  that  have  since  been  published  dealing  with  this  topic  can  be  found 
in  ABBASSI  (1964)  . 

The  essential  problem  of  significance  testing  is  the  selection  of 

an  effective  testing  procedure.  As  a  criterion  for  comparing  tests, 

NEYMAN  and  PEARSON  (1933)  introduced  the  notion  of  the  power  of  a  test. 

The  power  is  the  probability  that  the  test  will  reject  the  tested  hypothesis 

when  it  is  false.  Hence,  when  choosing  a  good  test,  we  seek  one  with  large 

power.  With  this  in  view,  the  object  of  this  thesis  is  to  derive  an 

approximate  expression  for  the  distribution  function  of  a  sample  serial 

correlation  coefficient  in  a  linear  Markov  process  and  thus  obtain  the 

approximate  power  of  a  test  of  serial  independence. 

Consider  the  linear  Markov  process  in  discrete  time  defined  by  the 

stochastic  difference  equation 

x  =  p  x  .  +  e 
s  s-1  s 

for  all  s  ,  where  the  real  constant  p  is  called  the  "lag  one  serial 

correlation  coefficient"  and  the  [  e  ]  are  independently  and  identically 

s 

distributed,  N(o  ,  l)  ,  random  variables.  In  the  present  thesis,  we  assume 

the  linear  Markov  process  to  be  stationary,  that  is,  I  P  ’l  <  1  .  Let 

x,  .....  x  be  a  sample  of  n  observations  from  the  process.  As  our 
1  n 

sample  estimate  of  the  serial  correlation  coefficient  p  we  take  the  ratio 


2 


(1.1) 


C 

A  O 

p  =  r 


where  and  C  are  the  following  quadratic  forms  in  the  observations, 


(1.2) 


and 


r  1  d 

C  =  2X1  +  X1X^  +  ...  +  x  -r  p 

o  1  12  n-1  n  ^ 


l  2  -2 

x  +  -p.  x  nx 

n 


2  2 

C  =  x  +  .  .  .  +  x 
1  n 


~2 

nx 


x  = 


X-  +  .  .  .  +  X 

1_ n 

n 


PATTON  (I96I)  derived  the  approximate  distribution  of  the  sample 

serial  correlation  coefficient  defined  above  using  the  method  of  DANIELS 

(1956)  .  The  approximate  renormalized  density  function  was  found  to  be 

N 

v2v2  -2 


where 


,  ,an  K(n,p)(l  -  p  )  ,  ,  ”p4  s 

h(p)  -  S  " -  n7!  (  1  +  0  (■  )  ) 


(1  -  2PD  +  p2) 


and 


K(n>p)  =h~ 


r  (Ltij 

V  2  ' 


r(|  +  1)  [(1  -  p2)  (n  -  1)  +  4  ] 


2 


N  =  n  -  1  + 


1  -  P2  * 

This  result  was  used  to  obtain  an  approximate  expression  for  the 
distribution  of  the  ratio  of  the  square  successive  difference  to  the  square 


difference,  n-1 


D  =  i=l 


(x.  -  -  x. V 
v  l+l  1' 


n 


Z  (xi  '  5)2 


i=l 

The  statistic  D  is  related  to  p  by  p  =  1  -  ^  D  and  for  large  n  is 
similar  to  the  ratio  of  the  mean  square  successive  difference  to  the 


variance, 


. 


-  3  - 


n- 1 


n 


i=l 


n 


(n-l)V  ( 


i=l 


The  distribution  of  the  statistic  rj  has  been  discussed  by  VON 
NEUMANN  (1941)  and  (1942)  and  WILLIAMS  (1941)  ,  for  the  case  where  the 


samples  are  taken  to  consist  of  independent  normal  random  variables 


x1  ,  .  .  .  ,  .  Patton  was  able  to  obtain  an  approximate  expression  for 


the  distribution  of  q  by  letting  p  equal  zero  in  the  approximation 
derived  for  the  distribution  of  D  . 

In  the  present  paper  we  take  the  approach  of  PATRICK  (1964)  to  obtain 
an  approximate  expression  for  the  distribution  function  of  the  sample  serial 


correlation  coefficient  p  .  This  approach  incorporates  the  methods  of 


GURLAND  (1948)  and  DANIELS  (1956)  with  the  asymptotic  techniques  developed 
by  ERDELYI  and  WYMAN  (I963)  and  (1963a)  for  deriving  asymptotic  expansions 
for  complex  integrals  in  terms  of  Confluent  Hypergeometric  functions. 

Briefly,  to  Gurland's  method  of  deriving  inversion  formulae  for  characteristic 
functions  is  applied  Daniels'  method  of  steepest  descent  to  obtain  a  complex 
integral  representing  the  distribution  function  of  the  sample  serial  correlation 
coefficient.  The  techniques  of  Erdelyi  and  Wyman  are  then  used  to  obtain  an 
asymptotic  expansion  for  the  complex  integral  in  terms  of  the  Confluent 
Hypergeometric  function. 


Our  review  of  Patrick's  method  omits  his  development  for  deriving  the 


asymptotic  expansion  of  the  integral 


[L  ] 

In  the  present  paper  we  simply  assume  his  results,  which  may  be  found  in 


(l  -  t^)S  dt  . 


detail,  in  PATRICK  (1964)  ,  Appendix  III  . 


-  4  - 


Chapter  II  is  a  review  of  Patrick’s  approach  for  the  case  of  a  serial 
correlation  coefficient  defined  in  a  circular  universe,  and  in  Chapter  III 
we  obtain  an  approximate  expression  for  the  distribution  function  of  the 
sample  serial  correlation  coefficient  defined  above.  In  Chapter  IV,  we 
find  the  approximate  power  of  a  test  of  serial  independence,  based  on  our 
sample  estimate  of  the  serial  correlation  coefficient,  for  different 
alternatives  p  and  sample  size  n  .  These  results  are  used  to  obtain  the 
power  curves  for  the  various  sample  sizes. 
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CHAPTER  II 

THE  APPROXIMATE  DISTRIBTUION  FUNCTION  OF  AN 
ESTIMATE  OF  THE  SERIAL  CORRELATION  COEFFICIENT 
IN  A  CIRCULAR  UNIVERSE 


This  chapter  serves  as  a  review  of  the  asymptotic  approach  adopted 

by  PATRICK  (1964)  to  obtain  approximate  expressions  for  distribution 

functions  of  sample  estimates  of  serial  correlation  coefficients  in  linear 

Markov  processes.  For  this  purpose,  we  consider  Patrick's  paper  with 

particular  reference  to  his  development  for  the  circular  case. 

We  define  the  linear  Markov  process  by 

x  =  p  x  ,  +  e 
s  s-1  s 

for  all  s  ,  where  the  real  constant  p  is  called  the  serial  correlation 

coefficient  and  the  [  e  ]  are  independently  and  identically  distributed, 

s 

N(o,l)  ,  random  variables.  The  process  is  called  circular  if 

x  =  p  x  ,  +  e 
s  s-1  s 

for  s  =  1  ,  .  .  .  ,  n  ,  where  x  =  x  and  p  <  1  . 

on1' 

Let  x,  ,  .  .  .  ,  x  be  a  sample  of  observations  from  the  circular 
1  n  r 

process.  They  have  a  joint  multivariate  normal  distribution  since 
e  ,  .  .  .  ,  e^  are  independent,  N  (o  ,  l)  ,  random  variables. 

This  distribution  is  given  by 
(2.1)  dF  (X]L  ,  .  .  .  ,  xn) 

(1  -  Pn)  r  1  r/,  2  w  2  2v 

=  n  exP  [  "  g  +  p  )(xi  +  •  •  •  +  xn)  " 

Yd 


(2it) 


2°  (xx  x2  +  x2  +  .  .  .  +  xn_t  xn  +  xn  xx))]  dxx 


.  dx  , 
n 


which  may  be  written  in  matrix  notation  as 


(2.2) 


dF 


(x  ,  .  .  .  ,  x  )  =  ( 1  exp  [-  \  x'  i~l  x] 

1  n  (27t)n/2  2 


dx,  .  .  .  dx  , 
1  n 


where. 


6  - 


xf  =  (xx 


X  ) 

n 


and 


i 


-i 


1+P2 

-p 

0  ... 

0 

0 

-p 

-p 

1+p2 

-p  ... 

0 

0 

o  1 

0 

-p 

1+p2  .  .  . 

0 

0 

0  i 

0 

0 

-p  ... 

0 

0 

0 

0 

0 

0  ... 

-p 

0 

0 

0 

0 

0  ... 

:u 

a 

+ 

r— 1 

-p 

0 

0 

0 

0  ... 

-P 

1+p2 

-p 

-p 

0 

0  ... 

0 

-p 

1+p2 

with 


t"l\?  =  1  -  p 


n 


As  an  estimate  of  p  ,  we  take  the  function 

(2.3)  p(xL  ,  • 


•  f 


X„>  - 


X1X2  +  X2X5  + 


+  X  X  +  X  X 

n-1  n  n  1 


X.  +  xr 


9  9  9 


+  X  ,  +  X 

n-1  n 


which  in  matrix  notation  is  given  by 

(2.U) 


A  X*  N  X 

P  ”  x'  D  x 


where y 


N  = 


and 


D  = 


0 

1 

2 

0 

0 

♦ 

9 

• 

0 

0 

0 

1 

2 

1 

2 

0 

1 

2 

0 

9 

9 

• 

0 

0 

0 

0 

0 

1 

2 

0 

1 

2 

o 

9 

e 

0 

0 

0 

0 

0 

0 

1 

2 

0 

• 

9 

a 

0 

0 

0 

0 

0 

• 

o 

9 

V 

o 

o 

9 

9 

9 

9 

0 

0 

0 

0 

9 

9 

C 

0 

1 

2 

0 

0 

0 

0 

0 

0 

o 

9 

9 

1 

2 

0 

1 

2 

0 

0 

0 

0 

0 

o 

e 

0 

0 

1 

2 

0 

1 

2 

1 

2 

0 

0 

0 

o 

• 
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0 

0 

1 

2 

0 

~T 

0 

0 

0 

9 

e 

o 

0 

0 

0 

0~ 

0 

1 

0 

0 

9 

0 

9 

0 

0 

0 

0 

0 

0 

1 

0 

9 

9 

O 

0 

0 

0 

0 

0 

0 

0 

1 

9 

9 

• 

0 

0 

0 

0 

• 

9 

9 

• 

9 

9 

• 

9 

9 

• 

« 

0 

0 

0 

0 

9 

• 

0 

1 

0 

0 

0 

0 

0 

0 

0 

« 

9 

9 

0 

1 

0 

0 

0 

0 

0 

0 
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9 
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0 

0 

1 

0 

0 

0 

0 

0 

• 

9 

O 

0 

0 

0 

1 

r 

”  7  " 


Thus  N  is  an  indefinite  symmetric  matrix  of  constant  coefficients  ,  D 
is  a  positive  definite  symmetric  matrix  of  constant  coefficients,  and 

I  p  I  <  i  . 

If  we  denote  the  distribution  function  of  p  by  H  (•)  ,  then  by 
definition 

(2.5)  H  (r)  =  Pr  [  p  <  r  ] 


=  P 


x>  N  x 


<  r 


x  *  D  x 


where  r  is  any  real  number  such  that  I  r  I  <  1  .  In  the  following  discussion, 
for  reasons  that  will  become  apparent,  r  is  restricted  so  that  I  r  I  <  1  . 

Thus , 


x  ‘  N  x  . 

-=s*  .  —  —  <  r 


x  '  D  x 


=  Pr  [(xf  N  x)  <  r  (x  ’  D  x)  ] 
=  P  [  r  (N  -  rD)x  <  0  ]  , 


Putting 


x '  (N  -  r  D  )  x  =  u  (x  ,  .  .  .  ,  x  ) 


and  denoting  the  distribution  function  of  u  by  G  ( ♦ )  ,  we  have 
(2.6)  H(r)  =  Pr  [  p  <  r  ]  =  [  u  <  0]  =  G  (o)  . 

The  matrix  (N  -  r  D)  is  indefinite  symmetric  with  variable  coefficients. 


The  characteristic  function  of  u(x.  ,  ,  ,  .  ,  x  )  is 

1  n ' 


(2.7) 


z-'  00  00 

cpu(0  =  /  •  .  .  /  exp  [i^u(x1  ,  .  .  .  ,  x  )]  dF(x1  ,  .  .  . ,  x  ] 

*^-00  {  J  “00 


where  t  is  a  complex  variable,  and  u(x.  ,  .  .  „  ,  x  )  and  dF(x,  ,  .  .  .  ,  x  ) 

1  n  1  n 

are  defined  as  before. 

We  determine  the  distribution  function  of  u  from  its  characteristic 
function  by  considering  the  inversion  formula 


G(v)  -\-+r  /  [-Ijv]  q.u(5)  df 


(2.8) 


-  8  - 


which,  for  v  =  0  ,  gives  the  distribution  function  of  p  ,  that  is, 


(2.9) 

where 

(2.10) 


H<r>  - 1  ■  SI  /  d?  ' 


\ 


/  =  u?€  +  X  ^ 

e->0 


the  path  of  integration  being  the  line  Im£  =  0  . 

In  matrix  notation  the  characteristic  function  is 


(2.U) 


<PU  (?) 


 lil 


1. 

2 


/oo 

#  #  # 

'00 


,  00 


-00 


(a* W 8 

X  exp  [-|  x'  t  ^  x]  dx^  .  .  .  dx^ 


exp  [ i£x  '  (N  -  rD )  x 


111"2 


(at) 


n/2 


>  00 


-00 


exp  [  x'  {  t 

J  -00 


dx1  .  .  .  dx 
1  n 


There  is  an  orthogonal  transformation  that  will  diagonalize  {  £  ^  -  2i £ 
Performing  this  transformation  and  integrating  we  get 


(2.12) 


where 


q>  (5)  = 


111 


l. 

’2 


U 


If1  -  2i^  (M  -  rD)l  2 


It”1  -  2i^  (N  -  rD)  |"2 


] 


2i£  (N  -  rD))x  ] 


(N  -  rD))  . 


is  the  Jacobian  of  the  transformation. 
By  equation  (2.14)  we  have 


9  - 


(2.13) 


N  -  rD 


-r  i- 
r  2 

i  -r 

0  i 


0  0 
0  0 
i  o 


0 

_1 

2 


0 

0 

0 


0 

0 

0 

• 

-r 

X 

2 


0 

0 

0 


l. 

2 


1. 

2 

0 

0 

• 

0 

1 

2 


0  \  -r 


Thus 

(£.11*) 


|i_1  -  £i?  (N  -  rD)| 


l+p2+2ir£ 

-(p+iO 

0  ... 

0 

0 

- (P+i? ) 

-(p+i£) 

2 

1+p  +2ir£ 

-(p+i£)  •  .  • 

0 

0 

0 

0 

"(P+iO 

l+p2+2ir£  .  .  . 

0 

0 

0 

0 

0 

0  ... 

l+p2+2ir£ 

-(p+i£) 

0 

0 

0 

0  ... 

-(p+iO 

l+p2+2ir£ 

-(p+i?) 

-(0+1?) 

0 

0  ... 

0 

-(P+iO 

2 

1+P  +2ir£ 

(p+i0n 

l+p2+2ir£ 

P+i£ 

-1 

0  ... 

0 

0 

-1 

-1 

l+p2+2ir£ 

P+i£ 

-1  ... 

0 

0 

0 

0 

-1 

l+p2+2ir£  .  .  . 

P+i£ 

0 

0 

0 

0 

0 

0  '  *  * 

2 

1+P  +2ir£ 

P+i£ 

-1 

0 

0 

0 

0  ... 

-1 

l+p2+2ir£ 

P+i£ 

-1 

-1 

0 

0  ... 

0 

-1 

l+p2+2ir£ 

P+i£ 

=  (p+iPlB]  • 

Thus,  equation  (2.12)  becomes 


n 


-  7~Ifyn/2  jip/2  • 


(2.15) 


10 


Evaluating  |  B_  |  as  a  circulant  we  have 


(2.16) 


Ul  = 


/ ,  n,2 

 ( 1  ~  z  ) 


n 


where, 

(2.17) 

Hence, 

(2.18) 

and 

(2.19) 

where, 

(2.19a) 

Now, 

(2.20) 


z  + 


1  1  +  p  +  2ir £ 

2  P  +  a 


(1  -  2pr  +  pd) 
P  +  1?  -  2  (1  -  2rz  l  «ti 


<!>,.(?(*))  = 


1  -  P 


n 


u 


(1  -  2pr  +  p  ) 


2sn/2  ■ 


£(z)  =  ip 


iz( 1  -  2pr  +  p  ) 

1  -  2rz  +  z2 


(l  -  2rz+  z  ) 

/  nv 

(1  -  z  ) 


2,n/2 


C ' (z)  =  (1  -  2pr  +  PP)(1  -  Z  2) _ 

£  (2)  (1  -  2rz  +  z2)(l  “  Pz)(z  _  p)  ’ 


thus,  from  equations  (2.19)  and  (2.20)  we  have 


(2.21) 


<Pu(5(z))  £  (z) 
£  (z) 


n 


_ 1  -  P 

/-  2xn/2  -1 

(1  -  2pr  +  P  )  / 


( 1  -  z  ) ( 1  -  2rz  +  z  ) 

(z  -  p)(l  -  P z ) ( 1  -  zn) 


2vn/2  -  1 


some 


Since  we  are  prepared  to  tolerate  errors  of  magnitude  0(An) 

A  |  <1  ,  for  neither  r  nor  p  near  +  1  ,  we  omit 

n 


X-~— n  in  (2.21)  . 


1  -  z 


for 


11 


Thus ,  we  have 


(2.22) 


<PU(5(«))  t'(z) 


TTzT 


/i  2\  / .  _  2\n/2  ~  1 

( 1  -  z  ) ( 1  -  2rz  +  z  )  / _ 

(1  -  2pr  +  p2)"/2  _1  (z  -  e)(l  -  pz) 


(8.23) 

where, 


Now,  equation  (2.22)  is  of  the  form 


<PU(5(*))  £'(*)<**  K,  p 

—  ~  / - \  g(z)(l  -  2rz  +  z)  dz 

(z  -  p)  ’ 


S  (*) 


1)  £  =  £(z)  represents  the  transformation  (2.17)  , 

2)  K'  is  a  function  of  (p,r), 

3)  g(z)  is  an  algebraic  function  taking  real  values  for  real  z 
and  is  regular  in  a  simply-connected  subset  of  |  z  |  <_  1  ,  which  contains 
z  =  p  and  z  =  r  , 


4)  s  is  a  real  variable  taking  large  values  [s  =  ^  -  1]  . 

In  equation  (2.23)  we  have,  for  our  case, 

n 


and 


K'  =  (1  -  2pr  +  02)"  2  +  1 


 (1  -  •*) 


*<*>  ■ 

The  transformation  (2.I7)  maps  the  £ -plane  cut  from  -i  00  to 
-i  (1  -  p^)/2(l  -  r)  and  from  i(l  +  p)^/2(l  +  r)  to  +i  00  onto  the 
open  disc  (  z  |<  1  ,  with  the  cut  being  mapped  onto  the  unit  circle  |  z  j 
If  we  denote  the  z-image  of  Im£  =0  by  P  ,  then  P  represents 
transformed  path  of  integration  in  the  z-plane,  with  the  end  points 


z  _  r  _  i ~ 2  and  z  =  r  +ij  ^  ^2  corresponding  to  the  points  £  = 

and  £  =  00  respectively.  The  z-image  of  the  point  £  =  0  is  the  point 


=  1  . 
the 

-  00 


z  =  p  at  which  the  path  P  crosses  the  real  axis.  The  curve  P  is  symmetric 
about  Im  z  =  0  ,  thus  points  below  Im  z  =  0  lying  on  P  are  complex 


conjugates  of  points  on  the  curve  above  the  line  Im  z  =  0  . 

Thus  the  z-correspondent  of  the  ^-integral  in  the  inversion  formula 
(2.9)  is  defined  to  be 


(2.24) 


K'  r 
2lT1  <_/  [r  ] 


g(z)(l  -  2rz  +  z?)Sdz 


where, 

(2.25) 


1 

2rti 


5 


The  integrals  defined  are  curvilinear  on  T  ,  with  a  and  its  complex 

-  -)f  — -)f 

conjugate  a  points  on  T  as  are  z  and  z  .  Thus  the  z-integral 
exists  in  the  sense  of  a  Cauchy  Principal  Value  as  defined  in  Appendix  I  . 
Now,  it  is  clear  that 

1  g(z)  ]  <  00  for  z  =  r  +  ij±_r2 

Thus 

(2.26)  lim  _  .  1  ni  g(z)(l  -  2rz  +  z2f  =  0  , 

+  VT7  '  '  p) 


and 

(2.27)  lim  _  1  n'  g(z)(l  -  2rz  +  zS)S  =  0 


Thus  we  have 

(2.28) 


4 

J(r] 


lim 
a — >p 


rA 

/  +  / 

J. 

yj  J 

a 


A  =  r  +  -  r2 


where 


-  13  " 


We  wish  to  evaluate  the  integral 
(2-29)  ,f[r]  irhy  g(z)(l  -  2rz  +  Z2)s  dz  , 


where, 


Jr 

'I  l r] 


and 


r  p  a  nA  1 

=  lim 

/ 

+ 

r 

J  - 

J  J 

a  — ^p 

A 

a 

is  regular  in 

z  |  <  1  . 

>  <A  -  r  +  i/i  .  r2  ) 


We  deform  the  path  of  integration  to  become  the  straight  line  running 
from  z  =  r  +  i 


L/l  - 


T  to 


z  =  r 


1  -  2rz  +  z2  I  <  (l  -  r2) 


<  1  ,  except  at  the  point 


iL  r2  ,  and  passing  through  the 

critical  point  at  z  =  r  .  The  deformed  path  is  made  in  such  a  manner  that 

(2.50) 

for  all  z  on  the  deformed  path,  which  we  denote  by  L  . 

The  integral  is  regular  for  all  z  in 
z  =  p  ,  and  therefore  the  integrand  has  only  one  residue  term.  This  residue 
term  being  the  one  corresponding  to  the  point  z  =  p  . 

We  have  three  cases  for  consideration.  They  are  r  <  p  ,  r  =  p  ,  and 
r  >  p  .  At  present  we  consider  the  two  cases  r  <  p  and  r  >  p  ,  dealing  with  the 
case  r  =  p  at  a  later  stage. 

We  now  indicate  the  following. 

(i)  r  is  the  deformed  path  joining  A  a  p  a  A  .  (as  before) 

(ii)  L  is  the  straight  line  joining  A  r  A  . 

(iii)  C  is  the  curve  a  a  (the  indentation  around  z  =  p]  « 


(iv) 


X. 


[  ] 


as  before,  represents  the  integral  from  x^  to  x^ 


on 


the  appropriate  curve  signified  in  [  ]  . 

For  r  <  p  ,  the  line  of  integration  is  in  the  positive  direction 
along  the  simple  closed  contour  A  a  a  A  r  A  ,  which  excludes  the.  point 


-  13a  - 


r  >  p 


Fig.  0 

Admissable  Deformed  Paths  of  Integration 


14  - 


z  =  p  by  the  indentation  C  .  Thus, 


a 


f 

J  A 

[r] 


J  a 

[C] 


J  a 

[r-] 


+ 


J  A 

[L] 


0  , 


a 


Va 

[r] 


+ 


J  a 

[r] 


7A 

[l] 


+ 


a 


a 

[C] 


Hence, 


lim 

r 

a — >  p 


a 


r 

J  A 

[r] 


a  J 


[H 


r 

Ja 

[L] 


+  lim 
C 

a — ^P 


Now, 


a 


lim 

C 

a — >  P 


and  therefore. 


a 


(7-  p)  s(z)(l  -  2rz  +  z2)S 


=  Jt 


2\S 


i  g(p)(l  -  2pr  +  p  )  , 


(2.31) 


i!r  /rrl  irryy  g  (z)(1  ‘ 2rz  +  z2)S  di 


[r] 

=  g(p) ( 1  -  2pr  +  p2)S 


+ 


~r  f  7—^ — T  g( z )  ( 1  -  2rz  +  z2)s 

2iti  J  r  (z  -  P)  /v  ' 


A 

[L] 


Now  for  r  >  p  ,  the  contour  encloses  the  pole  z  =  p  ,  and 
is  taken  in  the  negative  direction.  Thus 


a 


[H  [c] 


nA  n  A 

+  /  +  / 

'J  a 

[r]  [l] 


=  -2*ig(p)(l 


so  that 


[C] 


dz  . 


the  integral 


.  2vS 

2pr  +  p  )  f 


-  15  - 


so  that 


Thus 

(2.32) 


r>  3.  r>  A 

I  +  I 

'*  A  'J  a 

[H  [r] 


K'  £  1 

2*1  J  (■  r]  (z“p) 


r 

[L] 


J 


r 


a 


2X  s 


2itig(p )  ( 1  -  2pr  +  p  )  , 


a 

[C] 


g(z)( 1  -  2rz  +  z2)S  dz 


=  “f  g(p)(l  -  2pr  +  p2)S 


sfi  J-k  g(z)(1  • 2rz  +  z2)S  dz  • 


A 

[L] 


By  PATRICK  (1964)  ,  Appendix  II 
K'  g(p)(l  -  2pr  +  p2)S  =  1 

so  that,  for  r  <  p,  (2.3I)  becomes 

1 


(2-33)  ifi  J  (Z^T 

'  [r] 


g(z)(l  -  2rz  +  z2)S  dz 


1  K'  r  1 

2  +  2jti  ,J  (z-p) 

A 

[L] 


g(z)(l  - 


2rz  +  z2)S  dz 


and,  for  r  >  p  ,  (2.32)  becomes 


(2.34) 


k'  r 

2«1  J 

in 


1 

2 


(z-p)  §(z)(1  "  2rz  +  z2)S  dz 


4* 


r 


2\S 


2ici  Jl  (z-p) 
[L] 


g ( z ) ( 1  “  2rz  +  z  )  dz  „ 


Thus  the  inversion  formula  (2.9)  (using  (2.24))  for  r  <  p  becomes 


(2.35)  H(r)  = 


K* 

2*i 


g ( z )  ( 1  "  2rz  +  z2)S  dz  , 


i  A 
[L] 


(p-z) 


and,  for 


r  >  p 


becomes 


16 


(2.36)  H(r)  =  1  -  g(z)(l  -  2rz  +  z2)S  dz  , 

[L] 

where  L  is  the  straight  line  from  A  to  A  . 

The  integral 

(2-57)  J-A  g(z)(1  '  2rz  +  z2)S  dz 

[L] 

where,  L  is  the  line  from  z  =  A  =  r  -  i  J  ^  _  ^2  to  z  =  A  =  r  +  ~ *2 

passing  through  z  =  r  ,  has  a  critical  point  at  z  =  r  .  We  change  the 
variable  of  integration  so  that  in  terms  of  the  new  variable  of  integration 
the  critical  point  is  at  the  origin.  We  use  the  linear  transformation 


(2.38) 

z  =  r  +  1CA-  r2 

so  that, 

(2.39) 

dz  =  i/l  -  r2  dt  ’ 

(2.^0) 

(p-z)  =p-r-it/1_r2 

(2.U1) 

(1  -  2rz  +  z2)  =  (1  -  r2)(l  -  t2)  , 

(2.42) 

g(z)  =  g(r  +  it j  L  _  r2  )  =  x.(t)  , 

and  L  ,  the  t-image  of  L  ,  is  the  straight  line  from  t  =  -1  to  t  +  1 
passing  through  t  =  0  [i.e.,  a  segment  of  Imt  =  0  ]  , 

Thus,  the  t-integrals  corresponding  to  (2.35)  and  (2,3 6)  are, 


for  r  <  p 
(2.14-3)  H(r)  = 


K'fl  -  r2)s+a 

2jt 


(p  -  r  - 


X(t)(l  -  t2)3  dt 


and,  for  r  >  p 


-  IT  - 


(2.44)  H(r)  =  1  - 


K'(l  -  r2)S+s 

2jt 


r1 

'[lh 


(r  •  P  +  £Vl  -  r2) 


x(t)d  - 


t2)Sdt  . 


(2.^5) 


where, 


We,  therefore,  seek  an  asymptotic  representation  for  the  integral 

P  +  1 


jJsl 


-1  (p  -  r  ’  ^  1  .  r2> 


(1  -  t2)s  dt 


[L*] 


(i)  L  is  the  line  Imt  =  0  , 

(ii)  s  is  a  real  variable  assuming  large  values, 

(iii)  p  ,  r  are  real  variables,  p  fixed  such  that  j  p  |  <  1  ,  and 
r  varying  continuously  over  the  set  {  £  l  |  £  <1  and  £  ^  p)  , 

(iv)  x(t)  Is  an  analytic  function  of  t  regular  in  some  neighbourhood 


x(t) 


< 


CO 


for  all  t  in  [-1  ,  +  l]  ,  i.e.  on  the  path 


of  t  =  0  and 
of  integration. 

Assuming  the  results  obtained  by  PATRICK  (1964)  ,  Appendix  III, 
we  have 


(2.1*6) 


r 


+i 


jM. 


i-i  p  - r  - 

[L  ] 


(1  -  t2)S  dt 


M  m 


P  -  r 

1  -  r‘ 


,(m,l) 


T  (m+k+-|) 


(0  - 


m=0  k=0 


1  -  r 


2  r 


X  sl'U(m+k  +  -J;m+k  +  -|;- 

1  -  r 

M  m-1 


(P  -  r 


c*”‘1,2)  r(m  +  k  +  •£) 


m=l  k=0 


1  -  r 


">  m+k-o 

} 


X  sky^m+k  +  i»m+k  +  -|;  ~ 

1  -  r 


18  - 


+  0 


J  (p  -  r)2 
i  !  -  r2 


r(M  +  |)  u  (M  +  | 


s(p  -  r)2  \ 

1  -  r2  /  _ 


where,  U(  ;  ?  )  is  the  Confluent  Hypergeometr ic  function  in  the  notation 
of  SLATER  (i960)  ,  and  the  ,  C^m  are  given  by 

rC  K. 

00  M 

(2.47)  £  £  C*"'1*  sk  t2^  =  i  [  x(t)  +  x(-t)] 

m=0  k=0 

2  2\ 

X  exp  [st  +  s  log  (1  -  t  )  ] 


and, 

(2.48) 


00  m-1 


y  y  c(tn-l,2)  gk  t2(m+k-l) 


m= 1  k=0 


g'  '  [x(0  '  x("t)]  exP  [st2  +  s  lo§  (!  -  t2)] 


2J  1  -  r‘ 

Now,  by  equations  (2.44)  and  (2.45)  >  we  have,  for  r  <  p  , 


(2.49)  H(r)  = 


K* (l  -  r2)S 

2* 


M  m 


cf,:L)  r  (m  +  k  +  i) 


m=0  k=0 


X 


(P  °  rV 

1  -  r2 


m+k 


k  1 1 
s  U 


m+k+pt  m+k  + 


s(p 


r)2 

2 

r 


~  1  M  m-1 

K'(l  -  r2)s+?  V  V 


2* 


_  (ni"*  1,2)  r1/  1  1  \ 

/  /  1  (m  +  k  +  ■§■) 

m=l  k=0 


X 


Ip_ 

1 


if 


,  m+k— g- 


U  m  +  k  +  2  ^  m  +  k  +  2 


+  0 


(p  -  rV 

1  -  r2 


M+i 


r  (m  +  |  )  u 


2  .  m  i  2  .  sL°-“.  r)‘ 


M  +  “  |  M  + 


1  -  r 


and,  for  r  >  p  , 


(2.50) 


H(r )  *  1  - 


2ti 


2,tS.  M  m  /  .  v 

(m 9 1 )  p  , 

C  ',T  (m  +  k  +  f) 


m  f  0  k~=  0 


-  19  - 


x 


(r  -  P) 

1  -  r2 


2  -'I  m+k 


U  +  k  +  -g-  ;  m  +  k  +  -g-  J  S^r  ^ 


1  -  r 


K 1 ( 1  -  r  )S+2 

2* 


2xS+i  ^  ^  (m-1,2) 


m=l  k=0 


r,  .  M  f(r  -  p)2r+k-g  k 
2(m  +  k  +  2)  |  PI  ® 

1  1  -  p  J 


X  u(  In  +  k  +  \\  m  +  k  +  -g-; 


1.  s(r  ~  P)‘ 


1-r 

2  ' 


+  0 


"(r  -  P)2  |M+4 


r(M  +  |  ) 


x  U  (m  +  |  ;  M  +  |  ;  ats^Ei 
V  1  -  r 


Now, 


(2.51)  r(«r)  -{^4 

1  1  -  r 


2  ^  0--1 


u  ff  ;  cr  ; 


s(p  -  r)‘ 


00 


P-1 


1-r  /Jo  J  (p  -  r V 


exp  [-sv]  dv, 


+  v 


1  -  r 


and 


1.1.  S(P  -  r)‘ 
2 


(2.52)  U  !  2>2»  1  o*  j  —  exP 


1-r 


s(p  -o2  r 

1  -  r2  ^s(p-r)2 
1-r2 


v  2  exp  [-v]  dv  . 


Using  (2»5l)  for  m  +  k  >  0  and  (2*52)  for  m  +  k  =  0  ,  we  get,  in  terms 
of  the  above  integrals,  for  r  <  p  , 

2  -1  nOO  1 


(2.53)  H(r)  = 


K*(l  -  r2)S  .(0,1) 


2* 


C^}^'  T  (|)  exp 


s(p  -  r) 

2 

-  1  -  r 


2  exp  [-v]  dv 


s(p-r 

1-r2 


+ 


K'(l  -  r2)S4p  -  r) 

2jf 


M  m 


Q  (m»  ^  ) 

k 


m=l  k=0 


<00 


m+k-g 


L J °  J  (p  -  XY, 


exp  [-sv]  dv 


+  v, 


.  K* (1  “  r  ) 


2vs+^ 


M  m-1 


2jt 


(m-1,2)  k 
Ck  S 


+  0 


00 


m=l  k=0 


M+i 


J 


r 


00 


1  -  r 


m+k.  **  2 


°  J  (p  -  r)‘ 


exp  [-sv]  dv 


+  v 


1-r 


J  (P  ~  ll 


1-r 


exp  [-sv]  dv 


20  - 


and,  for  r  >  p  , 


2\  s 


(2.54)  H(r)  =  1  -  K r  )  cC0*1)  r  (i) 

2it  o 


X  exp 


s(r  "  P) 

1  -  r2 


2  -i 


,  00 


v  2  exp  [-v] 


s(r-p)' 

2 

1-r 


„  Kl_(l  -  r2)S~2  (r  -  p) 


M  m 


2it 


(m,l)  k 
k 


a-  -  •  s 


m=l  k=0 


X 


,  00 


m+k-i 


o  f(r  -  p)‘ 


exp  [-sv]  dv 


+  v 


1-r 


K*  ( 1  -  r^)S+2 

2rt 


2  s+-  —  n^-1 


p(m-l,2)  k 
s 


m=l  k=0 


X 


,  00 


m+k-g- 


'  0  M 
^  1 


exp  [-sv]  dv 


+  v 


+  0 


00 


'  °  f(l  Z  p)‘ 


exp  [-sv]  dv 


+  v 


1-r 


Thus,  as  r  -*■  p  ,  (2,55)  becomes 


(2.55)  H(p  -  0)  = 


K*(l  -  p2)S  C*  (°’1) 
o 


„  0  i  M  m-1 

Els 

m=l  k=0 


*(m-l,2)  r  (m  + 


m 


+  0  1  M+i 


dv 


wH?r 


-  21 


and  (2.54)  becomes 


(2.56)  H(p  +  0)  =  1  - 


K*  (1  -  p2)S  C*(0,1) 


„  O  .,1  M  m"1 

K  P2)S+2  y  y  c*(m-l,2) 

m=l  k=0 


r(m  +^k  ° 


m- 


+  0 


15F 


where, 

(2.57) 

K  =  K'  at  r  =  p  . 

(2.58) 

C*^  =  C  at 

k  k 

and  the  interchange  of  limits  and  integration  is  justified  since  the 
integrals  are  uniformly  convergent. 

Now,  by  (2.47)  , 

(2.59)  c*(0,1)  =  X*  (0) 

where,  by  equation  (2.42) 

x*  (0)  =  8(P) 


and 

(2.60)  x(fc)  =  X*  (fc)  at  r  =  P  • 


By  PATRICK  (1964),  Appendix  II, 


(2.61) 


K  (1 


2xS 


£1 


(O)  K 


* 


(equation  15) 

(1  -  p2)S  g(p)  _  1 

o  “  2  * 


Thus , 

(2.62)  H(p  -  0)  =  H(p  +  0)  =  H(p) 

1  M  m- 1 

1  K  (1  -  p2)S+2  V  V  *(m-l,2)  r(m  +  k  -  k) 

=  2  '  2jt  L  L  k  '  m-i 


m=l  k=0 


22 


+  0 


55¥ 


We  take  equation  (2.62)  to  be  the  asymptotic  value  of  H(r)  at  r  =  p 
Thus,  in  equations  (2.53)  and  (2.54)  for  M  =  0  ,  we  have  the 
following  asymptotic  expressions,  for  r  <  p 

(2.63)  H(r)  =  K'(lg~  r  ^  X  (0)  T(i) 


X  exp 


S(P  ~  O 

1  -  r2 


2  -i 


,  00 


v  2  exp  [-  v]  dv 


‘(p-r)‘ 

1-r2 


,00 


1. 

.2 


+  0 


{^4*  •} 


exp  [-sv]  dv 


and,  for  r  >  p 


2\S 


(2.61*)  H(r)  =  1  -  X  (0)  r  (4) 


X  exp 


s(r.  -  P.) 

1  -  r2 


2 


,  oo 


v  2  exp  [-v]  dv 


?(r-  ~  P.)‘ 

i  2 

1  -  r 


00 


+  0 


JL 

.2 


Jo  f(»  -  p)2 

1  l-r2 


exp  [~sv]  dv 


+  v 


which,  for  r  =  p  ,  give 

(2.65)  H(p)  =  i  +  0  (  -T 

X  s2 

which  would  be  an  approximation  to  the  first  term  in  (2.62)  .  Thus,  we 
take  as  the  asymptotic  representation  for  r  <  p  , 


-  23  - 


2  \  s 


(2.66)  H(r)  ^  r(i) 


exp 


s.(.P„-_r) 

1-r2 


2  n 


s(p-r)' 

1-r2 


v  2  exp  [-v]  dv 


and,  for  r  >  p  , 


2Ns 


(2.67)  H(r)  ^  1  -  T(i) 


exp 


s(r  -  p) 

1  -  r2 


2 


r 


00 


J 


;(r-p)‘ 

1-r2 


v  exp  [-v]  dv 


Equation  (2.42)  yields 

(2.68)  x(0)  =  g(r)  » 

Thus,  corresponding  to  equations  (2.35)  and  (2.36)  ,  we  have  for  r  <  p 


2\  s 


(2.69)  H(r)  „  g(r)  r  (4)  exp 


s(r  -  p) 


2  -1 


r 


1-r  J  s(r-p)‘ 
1-r2 


v  2  exp  [ - v ]  dv  , 


and,  for  r  >  p 


2vS 


(2.70)  H(r)  ^  1  -  Kt"^1  ~2k  r(i) 


exp 


s(r  -  P) 

1-r2 


2 


00 


f  0  v  2  exp  [-v]  dv 
J  s(r-p) 

1-r2 


Now,  by  equations  (2.22)  and  (2.23)  we  have 


(2.71) 


(i)  s  =  |  -  1 


(ii)  K'  = 


(1  -  2pr  +  p  J2 


?  -  1 


( iii)  g(z)  = 


1  -  z 
1  -  pz 


Thus  equations  (2.69)  and  (2.70)  give,  for  r  <  p 
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n 


(2.72)  H(r)  ^  ^ 


,,  2\2 
U  ~  r  ) 


P  n  -  1 

(1  -  2pr  +  p2)^ 


exp  L 


(1  -  1) 


2  -1 


1  -  r 


(1  -  pr) 


X 


(|  - 1)  ^^4  1 

1  -  r 


v  2  exp  [-v  ]  dv 


and,  for  r  >  p 

(2.73)  H(r)  ~  1  -  57^ 


n 


/ 1  2\2 

( 1  -  r) 


2  ^  —  X 

(l  -  2pr  +  p 


>  00 


/SL  _ 

K2  '  2 
exp  L  1  -  r 


2 


(1  -  Pr) 


X 


S-  1) 


(r  -  pV 

1  -  r2 


v  2  exp  [-v  ]  dv 


and,  in  particular,  for  r  =  p  , 


(2.74)  H(p)  \  . 


The  magnitude  of  error  incurred  by  the  above  approximation  is 


(2.75) 


R  =  0 


n 


as  n  — ►  00 


This  order  relation  being  uniform  in  r  and  n  , 
Now,  for  the  case  r  <  p  ,  consider 


(2.76)  exp 


Then, 


n  _  1  ^  (P  "  r) 


2 


1  -  r 


> 00 


v  2  exp  [-} 


n 

2 


-  1 


(p  -  r.r 

1  -  t‘ 


2 

Put  V  =  \  t 


dv  =  tdt  , 


v"2  =  v 2/t 


]  dv  . 
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and. 


v  = 


n 

2 


corresponds  to  t  =■ 


\l n-2  (p  -  r) 

Jl  -  r2 


Thus, 


(2.77 ) 


_i  /  Vn-2  (p  -  r) 


n/1  - 


\2  -i 


J 


$ 

-  \fn-2  (p  -  r) 

i 

• 

•-t 

ro 

T  _ 

9 

-  \l n-2  (p  -  r) 

■J  1  -  r2  . 

where  0  (  )  is  the  standardized  normal  distribution  function  and  cp  (  ) 

the  standard  normal  probability  density  function. 

The  corresponding  result  for  r  >  p  is 


(2.78) 


0 

-  s/n-2  (r  -  p) 

■Ji  -  . 

9 

 -  n/ n-2  (r  -  p) 

«/l  -  r2 

Thus,  by  equations  (2.72)  and  (2.73)  the  approximate  distribution 
function  for  p  is  given,  for  r  <  p  ,  by 


n 


H(r) 


M  2x2 

(1  "  r  ) 


<D 


v  n°2 


P  -  r 


Jl  r^” 


(1  -  pr )  ( 1  -  2pr  +  p2)2  1  9 


J 


1  -  r 


(2.79) 
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and,  for 


(2.80) 


r  >  p 


n 


H(r) 


1  > 


N^rt 


/ ,  2*2 
(1  ~  r  ) 


P  -  -1 

(1  -  pr ) ( 1  -  2pr  +  p  )2 
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CHAPTER  III 

THE  APPROXIMATE  DISTRIBUTION  FUNCTION 
OF  AN  ESTIMATE  OF  THE  SERIAL  CORRELATION 
COEFFICIENT  IN  A  STATIONARY  LINEAR 
MARKOV  PROCESS 


For  the  linear  Markov  process  defined  in  Chapter  II,  we  derive 
an  approximate  expression  for  the  distribution  function  of  a  sample 
estimate  of  the  serial  correlation  coefficient  p  with  unknown  mean 
in  the  stationary  case.  Our  sample  estimate  p  of  p  is  given  by 
(1*1)  and  (1.2)  ,  for  which  the  approximate  distribution  was  obtained 
by  PATTON  (I96I)  . 

We  assume  the  linear  Markov  process  to  be  stationary,  that  is, 


(3*1) 


x 


=  P  X, 


,  +  e 
■1  s 


for  s  =  1  ,  .  *  *  ,  n  ,  where  |  p  |  <  1  and 


cov  (x 

s 


o 


Let  x^  ,  .  .  *  ,  x^  be  a  sample  of  observations  from  the  stationary 
process.  They  have  a  joint  multivariate  normal  distribution  since  e  ,  ,  , 
e^  are  independent,  N(0  ,  l)  ,  random  variables. 

This  distribution  is  given  by 


dF  ( 


x  ) 

n' 


2  — 

^  '  ^2  exp  [Xj  +  (1  +  pO(x?  +  •••+xn) 

(2*)n/2  1  1  " 


+  x^  -  2p  (x-jXg  +  * 


.  +  x  .X  ) ]}  dx. . . ,dx  , 
n-1  n  1  n 


(3.2) 


y  •  •  ♦  9 
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which,  in  matrix  notation  may  be  written  as 


(3.3) 

dF 

(*]_» 

x )  _  (i  -  p¥ 
•"  * ‘ 

exp 

[- s  x'  r1  x 

]  dxx  .. 

where 

x' 

=  (xx 

»  •••  »  xn)  > 

and 

1 

-P 

0  ... 

0 

0 

0 

-P 

1  + 

2 

P  -P  ... 

0 

0 

0 

0 

-P 

.  2 

1  +  P  ... 

0 

0 

0 

0 

0 

-p  ... 

0 

0 

0 

ii 

1—1 

1 

■w-l 

0 

0 

0  v  «  • 

-p 

0 

0 

0 

0 

0  ... 

1  + 

P2  -p 

0 

0 

0 

0  ... 

-p 

1  +  p2 

-P 

0 

0 

0  ... 

0 

-p 

1 

and 

1  i 

1 

1 2 

=  (1  -  p2)2  . 

When 

the  mean  is 

unknown  we  take  the 

ratio 

(3.^) 

A 

P 

c 

c 

,  as  an  estimate  of 

P  > 

where  , 

.  dx  , 
n 


—  2  —  w 

C  =  i  (x1  -  x)  +  (x1  -  x)(x2  -  x)  +  . 


+  (x  -  x) ( X  -  x) 
n-1  '  n 


+  i  (xn  -  x)2 

-2 

=  c  -  nx  , 


Co  =  (k1  ‘  + 


.  +  (x  -  x)2 
'  n  ' 


-  c  -  nx 


(3.5) 
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with 


c  = 


1  2  i  2 

2  X1  +  XXX2  +  •••  +  Xn_lXn  +  2  Xn 


c  = 


X,  +  ...  +  X 
1  n 


and , 


x  = 


X  +  ...  +  X 

1_ n 

n 


In  matrix  notation 


(3.6) 


A  x'N  x 

x'D  x 


where. 


N 


(i  -  £>  (I  -  i) 


a 


b  <*  -  b 

(i  -  b  (4)  (4  -  b 


(-b 


/ 1  1  \ 

/  lx 

/ 1 

lx 

<5  -  n> 

™n^ 

(2 

'  n 

(i  - 

b 

4) 

n' 


(4  -  b  (i  -  b 


and 


9  |i— * 
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D 


with 


(1  -  i) 

n 


(1  “  ~) 
n' 


(1  -  ~) 
v  n 


(.i) 


(1  -  i) 

n' 


(i  -i) 


(i 


A  , 

P  |  <  1  . 


A  f  N 

Denoting  the  distribution  function  of  p  by  H  (*)  ,  the  following 
results  of  Chapter  II  are  again  valid. 

(3.7)  (i)  H  (r)  =  Pr  [p  <  r]  =  Pr  [u  <  0]  =  G  (o)  ,  where 

G  (*)  denotes  the  distribution  function  of  u  (xf  ,  ...  ,  x^)  =  x'  (N  -  rD)  x 


(ii)  The  characteristic  function  of  u  (x-  ,  ...  ,  x 


n 


)  is 


cp  (5)  =" 

u 


1_ 

■ 2 


If1  -  2i£  (N  -  rD)  |  2 


where 


.-1 


t  -  2i£  (N  -  rD)  |  '2 


is  the  Jacobian  of  the  transformation. 
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(iii)  The  distribution  function  of  p  is 

i  r  ^  (S) 

H  (r)  -  5  -  ,f  ~T~  d  l:  ’ 


where, 


£ 


,i 


=  lim 
T  — ►  00 

G  — ►  0 


r 


-e 


J  -r 


f 


j 


the  path  of  integration  being  the  line  Im£  =  0  * 
From  equation  (3 .6)  we  have 

i_r+“(r-i)  i+r(r_1) 


(3,8)  (N  -  rD )  = 


1  1 


^(r-1) 


(^(r-D) 


”r+~(r“l)  i+“(r-l) 

4+^(r"l)  -r+”(r-l)  i+^(r-l) 


(”(r-D) 


r+~‘( r_1 )  i+“(r=l) 


Thus 


-1 


(3.9)  £  -  2i£  (N  -  rD)  = 


e  +  a  b  +  a 


b  +  a  c  +  a  b  +  a 


(a) 


b  +  a  c  +  a  b  +  a 


(a) 


c  +  a  b  +  a 


b  +  a  e  +  a 


* 


B 
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where, 


e  =  1  -  i£  (1  -  2r ) 
21  t  (l  -  r) 


a  = 


(3-10)  / 


n 


b  =  -  (p  +  i£) 

V_  c  =  1  +  p2  +  2i£r 
Thus,  by  equation  (3.7)  >  (ii) 


(3.11) 


<!>  (5)  = 

u 


a  -  P2)^ 


*.x 


B 


To  evaluate  |B  I,  we  observe  that  1 B I  is  of  the  form  of  the  determinant  R  in 
1  —  1  •  '  — '  n 

Appendix  II  .  By  equation  (11.14)  ,  we  have 


(3.12) 


z  + 


1  =  1  +  p  +  2itr 
z  P  +  i£ 


1  +  p  -  2pr  +  2pr  +  2i£r 

P  +  i£ 


Thus , 
(3.13) 


1  +  p  -  2pr 

P  + 


+  2r  . 


P  + 


.£  =  (1  +  P2  -  2pr) 

(z2  -  2zr  +  l) 


From  equations  (II.IO)  ,  (il.ll)  and  (IX.15)  ,  we  get 


(3.14) 


R  = 
n 


1  -  a 


2  ( e  -  b  -  c )  -  n  ( c  +  2b  )  [ 


(c  +  2b)' 


X 


(eX  .  -  b2X  )  +  a(e,  -  c  -  b.)  I  2X  +  2(-l)n“1bn-1  1 
n-1  n-2'  (c  +  2b)2  1  n-1  ' 


where. 


X  = 


L±1 


n 


2\  n 


n  (1  "  *  )* 


z2n  -  Z2  -  (z2n  -  1)  Z  (-  |  ) 
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We  can  now  evaluate  the  determinant  |B  \  by  making  the  necessary 
substitutions  into  equation  (3,1k)  from  equation  (3.I0)  . 

Now,  using  equations  (3.I0)  and  (3.I3)  , 

(3.15) 

1  +  P< 1  -  2r) 


£  =  1  -  (1  -  2r) 

b  (P  +  i£) 


P  +  i£ 


+  2r 


(1  +  p  -  2pr)(z2  -  2zr  +  l)  +  (2r  -  l)(l  +  p2  -  2pr)  z 

2 

z  (1  +  p  -  2pr) 

2  2 

(pz  -  l)(z  -  p  +  2pr)  +  z  +  1  -  2prz 

2  * 
z  (1  +  p  -  2pr ) 


and  equation  (3*13)  gives  , 

(3.16)  (-b)n  =  (p  +  i5)n 

( 1  +  p2  -  2pr)n  zn 

(z2  -  2zr  +  l)n 

Hence,  equations  (3. 14)  ,  (3*15)  and  (3. 16)  yield. 


(3*17) 


(1  +  p2  -  2pr )n 

(1  -  z2)(z2  -  2zr  +  l)n 


2 

z 


^  f(pz  -  l)(z  -  P  +  2pr )  +  z2  +  1  -  2prz2} 

(1  +  P2  -  2pr ) 


(1  +  P2  -  2pr )n~ 1 

(1  -  z2)(z2  -  2zr  +  l)n 


X 


(z2n  -  z2)(l  +  p2  -  2pr )  -  ( z2n  -  l){(pz  -  l)(z  -  p  +  2pr ) 


2  o  2  1 

+  z  +  1  -  2prz  J 
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We  also  have  the  following  required  results 


(3. 18)  e-c-b=l-i£  (1-  2r)  -  ( 1  +  p2  +  2i£r)  +  p  +  i£ 


2 

=  P  -  P 


(3.19) 


Equations 


(3.20) 


Also , 
(3.21) 


c  +  2b  =  1  +  p2  +  2i£r  -  2  (p  +  i£  ) 

=  1  +  p2  -  2pr  +  2(r  -  l)(p  +  i£ ) 

(1  +  p2  -  2pr)(z2  -  2z  +  l) 

=  o 

(z  -  2zr  +  1) 

(3. 18)  and  (3.19)  give 

e  -  c  -  b  _  (p  -  p2)(z2  -  2zr  +  l) 

C  +  2b  (1  +  P2  -  2pr ) ( z  -  l)2 

a  =  2  in”1  t,  (1  -  r) 

=  2  n’1  {(r  -  1)  p  -  (r  -  l)(p  +  i£)} 

2  (r  -  l)  (p  (z2  -  2zr  +  l)  -  z(l  +  p2  -  2pr)) 

n  (z2  -  2zr  +  l) 

2  (r  -  l)(p  -  z ) ( 1  -  pz) 

n  (z2  -  2zr  +  l) 


and,  from  equations  (3*18)  and  (3. 19)  > 

2  (e  -c  -b)  -  n  (c+2b) 

(c  +  2b)2 

(z2  -  2zr  +  l)2  f2(p  -  p2)  -  n( 1  +  p2  -  2pr)(z  -  l)2(z2  -  2zr  +  1)  *1 

(1  +  P2  -  2pr )2  (1  -  z)h 


(3c22) 
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thus , 


(3-23) 


And, 

(3.24) 


I  2(e-c-b)  -  n(c  +  2b)  1 

(c  +  2b)2  J 

2(r  -  l)(p  -  z )  ( 1  -  pz)  _  _ 1 _ _ 

n(z2  -  2zr  +  l)  ( 1  +  p2  -  2pr)2  (l  -  z)4 


X  (z2  -  2zr  +  l)2[2(p  -  p2)  -  n( 1  +  p2  -  2pr)(l  -  z)2(z2  -  2zr  +  l)  1 ] 


=  ^ — k)(g-“.-z).(l  Pz)  S2( P  -  p2)(z2  -  2zr  +  1) 

n( 1  +  p  -  2pr )  (1  -  z)4  ( 

-  n( 1  +  p2  -  2pr ) ( 1  -  z)2  [ 


e=l-i£(l-2r) 

=  1  +  (1  -  2r )  p  -  (1  -  2r)(p  +  i£) 

(1  +  p  -  2pr)(z2  -  2zr  +  l)  -  (l+p2  -  2pr)  z(l  -  2r) 

(z2  -  2zr  +  l) 

2  2 

(pz  -•  1 )  ( z  -  p  +  2pr )  +  z  +  1  -  2prz 

(z2  -  2zr  +  l) 


After  substituting  equations  (3*13)  >  (3*17)  >  (3*20)  ,  (3.23) 
and  (3.24)  into  equation  (3.14)  ,  and  omitting  terms  relatively 
0(n  ^)  we  have 


(3.25) 


(l  +  p2  -  2pr )n~3 _ _ 

(1  -  Z2)(l  -  z)2(z2  -  2zr  +  l)n  1 


-  2pr)}2  . 


X  {(1  -  p)(l  -  Pz)(l  +  p  +  pz 


-  z 
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Therefore,  substituting  equation  (3.25)  into  equation  (3, 11)  , 


we  get 

(3.26) 


where, 

(3.27) 


<PU  (£  (z))  ~ 


(1  -  Pg) 


pi 


(1  +  P  ”  2pr ) 


.2x4 


T2 


n-1 


X 


(1  -  z  )2(1  -  z)(z  -  2zp  +  1) _ 

(1  -  p z ) ( 1  -  p)(l  +  p  +  pz  -  z  -  2pr) 


£  (z)  .  ip  -  e  3*) 

(z  -  2zr  +l) 


Now, 


(3.28) 


£ '  (z)  =  (1  +  P2  -  2pr)(l  -  z2) 

^  ( Z ^  (z2  -  2zr  +  1 ) ( z  -  p)(l  -  pz) 


thus , 


(3.29) 


cpu(£  (z))  V  (z) 
£  (z) 


n-3 

(z2  -  2zr  +  1)  2 

n 


(1  +  P2  -  2pr )  2 


X 


(1  -  z2)V2(l  -  z ) ( 1  -  p2)s 


(1  -  pz)  (z  -  p)(l  +  p  +  pz  -  z  -  2pr ) { 1  -  p) 


As  seen  in  Chapter  II  the  transformation 


1 

z  +  - 
z 


1  +  P  +  2ir£ 

P  +  i£ 


2 

maps  the  £  -  plane  cut  from  -  ioo  to  -  i(l  -  p)  /2(l  -  r)  and  from 

2 

i(l  +  p)  /2(l  +  r)  to  4-  ioo  onto  the  open  disc  \  z  |  <  1  ,  with  the 
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cut  being  mapped  onto  the  unit-circle  |  z  |  =  1  . 
Letting 


(1  -  p)(l  +  p2  -  2pr)  2 


g(«)  -  - -  «) - 

( 1  -  pr )  (l+p+pz-z-  2pr ) 


> 


and 


s 


9 


equation  (3.29)  is  of  the  form 


(3.30) 


cpu  (£  (z))  V  (z) 

5  (z) 


dz 


(z  ^  p)  8  (z)(l  -  2zr  +  z2)S  dz 


where,  as  before 

1)  £  (z)  represents  the  transformation 

z  ,  I  =  1„+_P2.  -+-21_rl 

+  z  P  +  i£  ’ 

2)  K'  is  a  function  of  (p  ,  r)  , 

3)  g  (z)  is  an  algebraic  function  taking  real  values  for  real 
z  and  is  regular  in  a  simply-connected  subset  of  |  z  j  <  1  ,  which 
contains  z  =  p  and  z  =  r  , 

4)  s  is  a  real  variable  taking  large  values. 


Thus,  by  equations  (2.24)  and  (2.25)  the  z-correspondent  of  the 
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^-integral  in  the  inversion  formula  (3.7)  >  (iii)  is  given  by 


(3.31) 


K' 


f 


2«i  J  (z  -  p) 

[H 


g  ( z ) ( 1  -  2rz  +  z2)S  dz  , 


where 


J 


[H 


and 


lim 

r 

a—>  p 

r 


* 


+  r 

j . 


z  ->  r+r 


.717 


U) 


s  = 


"  -  5 


(ii)  K’  = 


O  i 

(1  "  P  )2 


ml 

(1  -  p)(l  +  p2  -  2pr)  2 


(iii)  g  (z)  = 


(1  -  zg)?/2  n  -  z) 


(1  -  pz)  (1  +  p  +  pz  -  z  -  2pr) 


Thus,  by  equations  (2*35)  and  (2.36)  the  inversion  formula 
(iii)  becomes,  for  r  <  p  , 


(3.32)  H  (r)  = 


K 


7 » 


2jti 


h  Trh y  g<z>(1 

[L] 


2rz  +  z^)S  dz 


and,  for  r  >  p 

(3.33) 


K' 


H  (r)  =  1  ’  aa  J-  (T- . p)  8  (z>(1  -  2rz  +  Z")S 

[L] 


(3.7), 


dz  , 
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where  the  path  of  integration  is,  as  before,  the  straight  line  from 
A  to  A  . 

With  the  integrals  (3.32)  and  (3.33)  satisfying  the  conditions 
of  Chapter  II  ,  the  general  results  previously  obtained  apply.  Thus 
equations  (2.69)  ,  (2,70)  and  (3.3I)  give  in  the  present  case, 
for  r  <_  p  , 


(3-31*) 


H  (r) 


u  p  ;  V- 


sii 

2v^c  (1  -  pr)2(l  +  p2  -  2pr )  2 


X  exp 


n.  Z-JL 

2 


(  P  -  r) 

1  -  r2 


2  n 


v  2  exp  [  ■ 


/n~3\  (p-r )' 

'  2  '  2 

1-r 


and,  for  r  >  p 


(1  -  p2)"2  (l-rl) 


(305)  H  (r)  ^  1  - 


2,n/2 


n 


2  n/x  (l  -  pr )2( 1  +  P2  -  2pr ) 


X  exp 


n 


2l  ^  (r  -  p) 


2 


1-r 


>  00 


r  v  2  exp  ( 


J  (  \2. 

.,n-3\  r-p ) 

(  o  )^T” 


1-r 


n  -3  / p  r 

Now,  by  equations  (2.77)  and  (2.78)  »  with  v  = 

1-r 


v]  dv 


v]  dv 


we  have,  for  r  <  p  , 
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(30  6) 


exp 


(n-“  3\(P  -  r) 

v  2  '  ,  2 

1  -  r 


2  i 


,  00 


v  2  exp  [-v]  dv 


/n~3\ (P~r  Y 

K  2  '  ,  2 

1-r 


$ 


=  /2 


<P 


-  \ln  -  1  (p  -  r) 


'n  -  (p  -  r) 

n/  1  - 


and,  for  r  >  p 


(3.37) 


exp 


rn  -  3x(r  -  P) 

v  2  '  ,  2 

1-r 


2  n 


,  00 


v  exp  [ - v ]  dv 


/n~3\ (r-p 

v  2  '  ,  2 

1-r 


0  "  -  ^~n  -  3  (r  -  p) 

.^2  — — 

-  Vn  -  5  (r  -  p) 


co 


\/l  -  r2 


Thus,  by  equations  (3*34)  and  (3. 35)  »  the  approximate  distribution 
function  of  p  ,  for  r  <  p  ,  is  given  by, 


(3.38)  H  (r)  ^  ^ 


(1  -  P2)'^  (1  -  rg1n/2 


(1  -  pr)2(l  +  P2  -  2pr) 


ElS 

2 


0 


X 


jg.  ■  3-Jp  ill 


<p 


-  sTn  -  5  I  P  -  r  | 

\/l  -  r^ 


and,  for  r  >  p  ,  by 
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(3-39) 


H  (r) 


!  (1  -  s^Z  (1  -  r2)"/2 

•m  p-r 

(1  -  pr)2(l  +  p2  -  2pr) 


-  42  - 


CHAPTER  IV 


THE  APPROXIMATE  POWER  OF  THE 
TEST  FOR  SERIAL  INDEPENDENCE 


For  the  linear  Markov  process  given  by  (3.1)  consider  the  hypothesis 
of  serial  independence,  that  is  p  =  0  .  To  perform  a  test  of  this 
hypothesis  we  require  a  suitable  testing  procedure.  That  is  a  practical 
as  well  as  dependable  method  for  deciding  to  accept  or  to  reject  the 
hypothesis.  As  a  criterion  for  selecting  such  a  test,  or  testing 
procedure,  NEYMAN  and  PEARSON  (1933)  introduced  the  concept  of  the  power 
of  a  test.  Simply,  this  is  the  probability  of  rejecting  the  hypothesis 
when  it  is  false.  Hence  a  desirable  test  is  one  with  large  power.  In 
this  chapter  we  investigate  the  approximate  power  of  the  test  of  serial 
independence  based  on  the  sample  estimate  of  the  serial  correlation 
coefficient  given  by  (3.4)  and  the  approximate  distribution  function 
obtained  for  this  estimate  in  Chapter  III.  The  approximate  power  is 
obtained  for  selected  values  of  the  serial  correlation  coefficient  p 
and  the  sample  size  n  ,  from  which  the  approximate  power  curves  for 
the  different  values  of  n  were  sketched. 

Consider  the  null  hypothesis  and  the  alternative  hypothesis 


(4.1) 


H  :  P  =  0 
o 

H1  :  P  /  0 


and  the  testing  procedure 
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(4.2)  (i)  Accept  if  |p  |  =  |p  ^  ,  x^)  |  >  c 

(ii)  Accept  HQ  if  |p  |  =  |p  (x;L  ,  ...  ,  x^)  |  <  c 

where  the  random  variable  p  is  the  sample  estimate  of  the  serial 
correlation  coefficient  defined  by  equation  (3.4)  ,  and  the  real 
number  c  is  determined  such  that 

Pr  [  |p  |  >  c  |  p  =  0  ]  =  .05  . 

The  power  of  the  test  for  the  alternative  p^  is  given  by 

(4.3)  P(PX)  =  Pr  [  | p  |  >  c  |  p  =  p1  ]  , 

and  the  power  function  of  the  test  is  the  function  of  p  ,  P(p)  . 

The  approximate  power  function  of  the  test  for  serial  independence, 

A 

in  terms  of  the  approximate  distribution  function  of  p  obtained  in 
Chapter  III  and  given  by  equations  (3.35)  and  (3.36)  ,  is  ,  for 
-c  <  p  <  c  , 


(b.k) 


p(p)  = 


■fs 


!jt 


fi  -  p2r^  a  -  c2)n/g 

Ell 

(1  -  pc)2(l  +  P2  -  2pc)  2 


0 

-  n/  n  -  3  1  P  -  C  1 

L 

-  sTn  -  3  p  -  c 

\/l  -  c^ 

(1  +  pc ) 2 ( 1  +  P2  +  2pc )  2 


X 
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$> 


X 


-  sfn  -  3 


P  +  c  1 


n/~1  - 


<P 


iJks-5- 1. 


P  +  c 


^yl_^ 


and  ,  for  -  c  <  c  <  Ip 


2n-4 


(4.5) 


p(p)  =  1  + 


(i  -  Ip  1  )  a  fi  -  c  ) 


.  ,V/2 


■JZn  Ell 

(1  -  c|p  |  )2(  1  +  |p|2+2c|p|)2 


X 


0 

cp 


-  N^n  -  3  |  |  p  1  +  c  | 


_  (l  -  Ip  1 2r*  (l  -  c2)n/2 

Ell 

(1  -  c  |p  | )2(1  +  |p  |2  -  2c  |p  | )  2 


<D 

"  ^~n  -  3  1  1  P  1  -c[ 

■ITT? 

<p 

-  \ln  -  3  p 

-  c 

where  c  is  such  that 


0.025 
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Table  1.  The  Approximate  Power  of  the  Test  for  Serial  Independence  at 
a  Significance  Level  of  5  $  . 


X 

10 

15 

20 

25 

35 

50 

0.0 

0.050 

0.050 

0.050 

0.050 

0.050 

0.050 

0.1 

0.057 

0.063 

0.069 

0.075 

O.087 

0.150 

0.2 

0.082 

0.106 

0.132 

0.157 

0.210 

0.288 

0.3 

0.129 

0.191 

0.253 

0.316 

0.434 

0.567 

0.4 

0.212 

O.334 

0.451 

0.530 

0.645 

O.790 

0.5 

0.351 

0.515 

0.589 

0.672 

0.811 

O.927 

0.6 

0.502 

0.587 

0.703 

O.78O 

0.918 

0.981 

0.7 

0.498 

0.660 

0.800 

0.890 

0.970 

0.996 

0.8 

0.473 

O.722 

O.87O 

0.943 

0.900 

0.999 

In  Table  1  we  have  tabulated  the  approximate  power  for  the  test 
(4.2)  ,  using  equations  (4.4)  and  (4.5)  >  for  selected  values  of  the 
serial  correlation  coefficient  p  and  sample  size  n  .  The  approximate 
power  function  is  symmetric  about  p  =  0  ,  and  therefore  we  have  calculated 
the  approximate  power  only  for  positive  values  of  p  . 

An  examination  of  Table  1  discloses  that,  for  a  sample  of  size  10  , 
the  approximate  power  decreases  for  values  of  p  greater  than  0.6  , 
instead  of  continuing  to  increase  as  would  be  expected.  Although  we 
have  not  been  able  to  pinpoint  the  exact  reason  for  this  we  suggest 
that  a  full  investigation  would  reveal  the  approximate  distribution 
function  for  p  ,  which  has  been  obtained  on  the  assumption  of  moderate 
p  and  r  ,  and  large  n  ,  suffers  a  considerable  loss  in  accuracy  here. 


'  M  .  J  t> 
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We  also  note  that  the  approximate  power  curve  for  n  =  15  ,  illustrated 
in  Appendix  III  ,  Fig,  1  ,  has  a  distinct  flattening  of  the  tails, 
whereas  the  curves  for  the  larger  values  of  n  do  not  exhibit  this 
characteristic.  This  would  indicate,  that  in  the  present  context,  a 
reasonable  test  could  not  be  based  on  a  sample  of  less  than  20 
observations . 

In  practise  a  test  is  usually  made  on  a  relatively  small  number  of 
observations,  generally  20  or  less.  This  is  due  to  the  time  and 
expense  that  would  be  otherwise  involved.  Now,  Table  1  indicates  that 
for  a  sample  of  20  observations  we  would  only  be  able  to  detect  a 
p  value  of  0,5  about  60$  of  the  time,  whereas,  for  samples  o  f  35 
and  50  observations  we  could  detect  this  value  approximately  80$> 
and  93  1°  of  the  time  respectively.  Even  a  p  value  of  0,8  would 
remain  undetected  of  the  time  while  a  sample  of  50  observations 

would  fail  to  detect  this  value  in  only  0.1 $  of  the  trails.  Hence, 
we  conclude,  our  testing  procedure  would  only  be  of  any  real  practical 
use  in  situations  where  large  samples  of  observations  could  be 


conveniently  applied. 
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APPENDIX  I 

CAUCHY’S  PRINCIPAL  VALUE  OF  AN  INTEGRAL 


Assumptions : 

(1)  The  function  f(z)  is  defined  in  some  region  of  the  z-plane. 

(2)  For  a  fixed  point  z  =  p  and  some  fixed  R  >  0  ,  f(z)  is 
regular  for  all  points  of  the  neighbourhood  |  z  -  p  |  <  R 
except  at  z  =  p  . 

(3)  The  point  z  =  p  is  a  simple  pole  of  f(z)  . 

(b)  C  is  a  continuous  curve  in  the  z-plane  passing  through  z  =  p 
and  intersecting  with  |  z  -  p  |  =  R  at  z  =  A  and  z  =  B  . 

Consider  the  integral 


(1.1) 


lim 
a  -*  p 
b  — ►  p 


f(z)  dz 


where  the  limits  a  ->  p  ,  b  -> p  are  taken  on  the  set  C  . 

We  wish  to  determine  sufficient  conditons  for  which  the  above 
integral  exists  in  the  sense  of  a  Cauchy  Principal  Value  . 

By  assumption 

f(z)  =  +  g(z) 

for  all  z  in  |  z  -  p  )  <  R  where  g(z)  is  regular  for  all  z 
in  | z  "  P  |  <  R  ,  and  K  is  constant  with  respect  to  z  . 


50 


Thus , 


r  n a 

rB 

- 

lim 

+ 

_ <1 

_ i 

/ 

a  -►  p 

'b 

b  — v  p 


f(z)  dz 


=  lim 
a  ->  p 
b  -»  p 


dz 


=  lim 
a  ->  p 
b  p 


dz  +  lim 
a  ->  p 
b  — >  p 


g(z)  dz  . 


Now,  g(z)  is  regular  for  all  z 


b  -*■  p 


in  |  z 


g(z)  dz 


P  |  <  R 


hence  , 

dz 


Also , 


lim 


a  -v  p 
b  ^p 


=  K  lim 
a  — >  p 
b  ->  P 


=  K  lim  log 
a  — ►  p 
b  p 


} 


K  lim 
a  — »■  p 
b  ->  p 
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=  K  log 


{ 


B  -  P 

A  -  p 


} 


b  — ►  p 


Thus , 


b  — +  p 


=  lim  K  log 
a  -►  p 
b  — v  p 


f(z)  dz 


+  K  log 


+ 


dz  . 


Hence,  a  sufficient  condition  that  (i.l)  exists  as  a  Cauchy  Principal 
Value  is  that 


(1.2)  lim  K  log  -<  V  is  finite, 

a  -►  p  J 

b  -►  p 


Consider  the  normal  to  the  curve  C  at  z  =  p  ,  and  the  transformation 


z  -  p  = 


te 


iO 


where  z  ,  t  and  p  are  complex  variables,  and  0  is  a  fixed  angle 
determined  so  that  the  normal  to  the  curve  C  at  z  =  p  in  the  z-plane 
is  transformed  to  be  the  line  Imt  =0  in  the  t-plane.  Thus  in  the  t-plane 
(1.2)  becomes 


lim  K  log  ~— 

t  ->  0  b 

a 


0 


A 
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where  t  ,  t  are  the  t-images  of  a  and  b  respectively,  and  the 

ci  D 

limits  are  through  the  set  C  (the  t-image  of  C  in  the  z-plane)» 

Now  consider  the  special  case  where  the  curve  C  is  symmetric  about 
Imt  =  0  .  In  this  case  t&  =  t^  (the  complex  conjugate  of  t^  )  , 


Thus  if  tb  =  |  tb  |  elCp 

0  <  cp  <  2* 


and 


lim 

t  0 
a 

t,  ->0 
b 


becomes  lim 

Cp  ->  0 


Thus  we  have 


lim  K  log 

t  —v  0 
a 

t,  — ►  0 
b 


lim  K 
cp  ->  $ 


=:  lim  K  log  (e"2lCp) 

cp  -►  0 

=  lim  -  2iKcp 

cp  $ 

=  -  2iK3> 


which  is  constant  and  finite  . 

[cf.  SANS ONE  and  GERRETSEN  (i960),  pp .  127  -  128] 
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APPENDIX  II 

EVALUATION  OF  A  DETERMINANT 


We  wish  to  evaluate  the  determinant 


(II. I)  R 


n 


e  +  a 
b  +  a 
a 
a 


(a) 


b  +  a 
c  +  a 
b  +  a 
a 


a 

b  +  a 
c  +  a 
b  +  a 


(a) 


b  +  a 
c  +  a 
b  +  a 
a 


a 

b  +  a 
c  +  a 
b  +  a 


a 

a 

b  +  a 
e  +  a 


This  may  be  written  as  the  (n  +  l)  X  (n  +  l)  determinant 


n 


(II. 2)  R 


n+1 


1 

0 

0 

0 


0 

0 

0 

0 


1 

e  +  a 
b  +  a 
a 


a 

a 

a 

a 


1 

b  +  a 
c  +  a 
b  +  a 


a 

a 

a 

a 


1 

a 

a 

a 


b  +  a 
c  +  a 
b  +  a 
a 


1 

a 

a 

a 


a 

b  +  a 
c  +  a 
b  +  a 


1 

a 

a 

a 


a 

a 

b  +  a 
e  +  a 


n+1 


If  we  subtract 


"a"  times  the  elements  of  the  first  row  from  each 


of  the  remaining  rows  we  get 
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i 

l 


J 

(n.3)  Rn+1  -  | 


+ 


1 

1 

1 

♦ 

♦ 

1 

1 

1 

-a 

e 

b 

* 

-• 

* 

0 

0 

0 

-a 

b 

c 

•tr 

0 

0 

0 

-a 

0 

b 

♦ 

♦ 

*• 

0 

0 

0 

-a 

0 

0 

• 

♦ 

♦ 

b 

0 

0 

-a 

0 

0 

V 

♦ 

• 

c 

b 

0 

-a 

0 

0 

-* 

* 

* 

b 

c 

b 

-a 

0 

0 

# 

* 

* 

0 

b 

e 

0 

1 

1 

•  *  • 

1 

1 

1 

-a 

e 

b 

•  •  • 

0 

0 

0 

-a 

b 

c 

♦  *  % 

0 

0 

0 

-a 

0 

b 

♦  ♦  ♦ 

0 

0 

0 

-a 

0 

0 

b 

0 

0 

-a 

0 

0 

♦  ♦  # 

c 

b 

0 

-a 

0 

0 

*  ♦  t 

b 

c 

b 

-a 

0 

0 

»  •  « 

0 

b 

e 

1 

1 

1 

»  •  ♦ 

1 

1 

1 

0 

e 

b 

»  t»  * 

0 

0 

0 

0 

b 

c 

0 

0 

0 

0 

0 

b 

♦  *  * 

0 

0 

0 

0 

0 

0 

•  •  « 

b 

0 

0 

0 

0 

0 

»  • 

c 

b 

0 

0 

0 

0 

♦  •  • 

b 

c 

b 

0 

0 

0 

♦  *  ¥ 

0 

b 

e 

n  +  1 


1 


a  +  1 
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or 


(II. 4) 


R  =  -  a  P  +  Q 
n  n+1  n 


where 


n+1 


and 


n 


0 

1 

1 

•  9 

• 

1 

1 

1 

1 

e 

b 

9  9 

9 

0 

0 

0 

1 

b 

c 

9  9 

9 

0 

0 

0 

1 

0 

0 

9  9 

0 

c 

b 

0 

1 

0 

0 

9  9 

0 

b 

c 

b 

1 

0 

0 

a  • 

0 

0 

b 

e 

e 

b 

0 

9  9 

9 

0 

0 

0 

b 

c 

b 

9  0 

♦ 

0 

0 

0 

0 

b 

c 

9  0 

e 

0 

0 

0 

0 

0 

0 

©  © 

© 

c 

b 

0 

0 

0 

0 

o  o 

• 

b 

c 

b 

0 

0 

0 

9  9 

• 

0 

b 

e 

sum  of 

rows 

2  to 

(n  + 

1) 

of  P  . 
n+1 

by 

(-!)/( 

n+1 


n 


and  adding  to  the  first  row,  we  have 
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-n 

1  e  +  b 

0 

0 

0 

l  e  + 

b 

c  +  2b 

c  +  2b 

<  «  « 

c  + 

2b 

| 

1 

e 

b 

♦  •  9 

0 

0 

0 

1 

b 

c 

•  •  * 

0 

0 

0 

1 

0 

0 

•  •  • 

c 

b 

0 

1 

0 

0 

•  •  •# 

b 

c 

b 

1 

0 

0 

•  *  t 

0 

b 

e 

Multiplying  the  sum  of  columns  2 

to 

(n  +  1) 

of 

Pn+1 

by 

(-l)/(c  + 

ab )  and  adding  the 

result  to 

the 

first  column 

,  we 

have 

2(e-b-c)  “  n(c+2b) 

,  e+b 

0 

0 

0 

,  e+b 

(c  +  2b)2 

c+2b 

•  •  * 

c+2b 

1  e+b 

c+2b 

e 

b 

•  •  • 

0 

0 

0 

0 

b 

c 

i  «  • 

0 

0 

0 

Pn+1  = 

0 

0 

0 

»  »  • 

c 

b 

0 

0 

0 

0 

«  »  • 

b 

c 

b 

■j  e+b 

c+2b 

0 

0 

» 

0 

b 

e 

2  ( e  -b  -  c)  -  n  ( c  +  2b ) 

(c  +  2b)2 


+ 


e  -  b  -  c 
(c  +  2b) 


M  + 
n 


e  -  b  -  c 
(c  +  2b) 


N 

n 


where. 
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(II. 6)  M  = 


n 


and 


(II. 7)  N  = 


n 


c+b-e 

c+2b 

0 

0 

•  •  9- 

0 

0 

c+b-e 

c+2b 

b 

c 

b 

9-9  9 

0 

0 

0 

0 

b 

c 

9  9  » 

0 

0 

0 

0 

0 

0 

c 

b 

0 

0 

0 

0 

9  9  9 

b 

c 

b 

0 

0 

0 

9  9  9 

0 

b 

e 

c+b-e 

c+2b 

0 

0 

9  9  9 

0 

0 

c+b-e 

c+2b 

e 

b 

0 

9  9  9 

0 

0 

0 

b 

c 

b 

9  9  9 

0 

0 

0 

0 

0 

0 

9  9  9 

b 

0 

0 

0 

0 

0 

9  9  9 

c 

b 

0 

0 

0 

0 

9  9  9 

b 

c 

b 

n 


n 


Expanding  in  terms  of  the  elements  of  the  first  row,  we  get 


(II. 8) 


M  =  c  + b e  x  +  (-i)"-1  c-+  > e  b"-1 

n-1 


n 


c  +  2b 


c  +  2b 


and 


(II. 9) 


N  .  c  +  b  :he  b"-1  +  (-i)"-1  c  +  b  -  e  x  , 

n  c  +  2b  c  +  2b  n-1 
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where, 


e 

b 

0 

♦  *  • 

0 

0 

0 

b 

c 

b 

•  *  ♦ 

0 

0 

0 

0 

b 

c 

•  ♦ 

0 

0 

0 

0 

0 

0 

•  •  ♦ 

c 

b 

0 

0 

0 

0 

♦  •  * 

b 

c 

b 

0 

0 

0 

•  *  » 

0 

b 

c 

Hence 


(II. 10) 


f  2(e  -  b  -  c)  -  n(c  +  2b)  "1 

^  (c  +  2b)2 


+ 


e  -  c  -  b)2 

(c  +  2b)2 


{ 


1 


+ 


2(-l)n'1  b"'1 


where , 

(II. 11)  Q  =  e  X  .  -  b2  X 

'  n  n-I  n-2 


Now,  from  equation  (II.9)  > 


(n.12) 


e 

b 

0 

•  •  ♦ 

0 

0 

0 

b 

c 

b 

•  *  ♦ 

0 

0 

0 

0 

b 

c 

♦  ♦ 

0 

0 

0 

0 

0 

0 

•  •  ♦ 

c 

b 

0 

0 

0 

0 

•  »  • 

b 

c 

b 

0 

0 

0 

0 

b 

c 
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e 

"  b 

-  l 

0 

♦  ♦  • 

0 

0 

0 

-  1 

c 

"  b 

-  l 

c 

•  •  • 

0 

0 

0 

0 

-  1 

"  b 

•  »  ♦ 

0 

0 

0 

0 

0 

0 

*  *  ♦ 

c 

”  b 

-  1 

0 

0 

0 

0 

•  9  •• 

-  1 

c 

b 

-  1 

c 

0  0  0  *  .  •  0-1 
=(-b)*An  . 


Thus,  if  is  the  determinant  consisting  of  the  first  j 

and  j  columns  of  A^  ( j  <  n)  ,  we  have 


(i)  Ax  =  - 


(ii) 


A„  = 


(§>(?>  -  1  • 


and  in  general 


(iii) 


A.  = 


"  b^  Aj-1  "  Aj -2 


The  solution  of  this  difference  equation  is 


(11.13) 


where, 

(I!. Ut)  z  +  I  =  -  t 


n 


rows 
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Hence,  equation  (II. 12)  gives 


(11.15)  xn  =  (-  b)n 


*  +  (g)  » 

1  -  z2 


}■" 


(-  b  ~  z)  2 

1  -  Z2 


-n 


b)n 


/,  2n  n 
(1  -  z  )  z 


>  .  22  .  (z2"  -  1)  z  (-  f) 


Thus,  from  equations  (il.lo)  ,  (il.ll)  and  (11.15)  ,  we  get 


(11.16) 


R  = 
n 


1  -  a 


2(e-b-c)-n(c+2b) 

(c  +  2b)2 


X 


e  (-  b) 


n-1 


/.  2\  11-1 

(1  -  z  )  z 


z2"-2  .  22  .  (z2n-2  .  1}  z  (.  1  ) 


(-I)"'1  bn  _ 

/ ,  2v  n-2 
(1  -  z  )  z 


z2n'4  -  z2  -  (z2"-4  -  1)  z  (-  §  ) 


2a(e-c-b)  ( -b ) 


(c  +  2b)' 


n-1 


/-  2\  n-1 

L(1  -  z  )z 


2n-2  2  /  2n-2  . n 

z  -  z  -  (z  -  1) 


X  z  (-  s  )  I-  +  (-I)"-1  b"-1 


. 
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APPENDIX  III 

POWER  CURVES 


This  appendix  serves  to  illustrate  the  data  of  Table  1  .  We 
reproduce  the  power  curve  for  the  cases,  n  =  15,  20,  25>  35  and 
50  ,  omitting  the  case  of  n  =  10  as  our  results  are  not  valid  here. 
The  curves  are  given  for  -0.8  <  P  <  0.8  ,  and  are  symmetric  about 
the  power  axis,  that  is,  the  alternative  p  has  the  same  power  as 
the  alternative  -p  ,  for  each  p  .  The  curves  give  us  an  overall 
indication  of  the  increase  in  power  as  the  serial  correlation 
coefficient  p  and  sample  size  n  become  larger.  Giving  us  an 
idea  of  the  effectiveness  of  the  test  for  various  n  . 


■ 


P  ( p  )Ho 


62 


o 


uO 

o 


4- 

'O 

o 


4 


■xj' 

o  o 


CN 

o 


•^T 

o 


-o 

o 


co 

o 


Figure  1  Power  Function  for  n  =  15 
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Figure  2  Power  Function  for  n  =  20 
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Figure  3  Power  Function  for  n  -  25 
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Figure  4  Power  Function  for 
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Figure  5  Power  Function  for 


